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Abstract 



The spin-orbit scattering effect on critical current through superconduc- 
tor/ferromagnet (SF) bilayers separated by an insulator (SFIFS tunnel 
junction) have been investigated for the case of absence of the supercon- 
ducting order parameter oscillations (thin F layers) . The analysis is based on 
a microscopic theory for proximity coupled SF bilayer for different bilayer 
parameters (boundary transparency, proximity effect strength, relative orien- 
tation and value of the F layers exchange field) . We find that the spin-orbit 
scattering considerably modifies dc Josephson current in SFIFS tunnel junc- 
tion. In contrast to a simple physical picture, the reduction of the exchange 
field effects is nonlinear in character getting its maximum in the field's region 
where the critical current enhancement or the transition to the 7r-state take 
place. Hence, for understanding the various experimental results on tunnel 
structures with thin F layers, the coupled effects of the exchange interaction 
and spin-orbit scattering must be considered. 

I. INTRODUCTION 

Nowadays progress in nanotechnology made it possible to produce nanostructures with 
new physical phenomena. This has led to renewed attention to hybrid systems consisting of 
superconductor (S) and ferromagnetic (F) metals, and displaying rich and elegant physics 
and having potential applications. Transport properties of SF structures with artificial 
geometry have turned out to be quite unusual. These have been treated by several authors 
[1-9], and the obtained results show that in ferromagnet the Cooper pair potential not 
only exponentially decays but also has an oscillatory character; i.e., in a ferromagnet the 
density of Cooper pairs is spatially inhomogeneous and the superconducting order parameter 
contains nodes where the phase changes by n. This causes the exchange field dependence of 
the Josephson coupling energy, and if the exchange energy H exc in magnetic layer exceeds 
a certain value, a crossover from 0-phase to 7r-phase superconductivity takes place. The 
phenomenon has been theoretically described for SFS weak links with thick F layer [1-5]. 
The crossover to the 7r-state even in the absence of the order parameter oscillations in thin F 
layer was also predicted for SFIFS tunnel junctions (where / is an insulator) with parallel 
alignment of layer's magnetization [6] . For an antiparallel orientation, the possibility of the 
critical current enhancement by the exchange interaction in SFIFS junctions with thin F 
layers was discussed both for small [7] and bulk [8] S layer thicknesses. Experimentally, 
the 7r-phase superconductivity has been observed for SFS weak links [9] and SIFS [10] 
and SNFNS [11] tunnel junctions with thick F layers, while the enhancement of the dc 
Josephson current has not been detected until now. 

Some features should be taken into account for theory to be adjustable for experimental 
results on SF structures, and the spin-flip processes are among the important ones. These 
processes can be induced, e.g., by the spin-orbit scattering centers presented in the film; 
other important source of the spin-flip processes for nanoscale hybrid structures is a strong 
electric field arising near metal-metal boundaries [12]. The basic physics behind spatial 
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oscillations of induced superconductivity in SF sandwiches with spin-orbit scattering has 
been recently discussed in [3,13]. As is known, in the presence of spin-orbit scattering the 
electron spin is no longer a good quantum number and the electron will change its spin state 
during a characteristic time Tso , while the Cooper pair will mix with its spin-exchanged 
counterpartner. This causes a pair to "feel" an exchange field which changes sign at a rate 
proportional to 1/tso , decreasing the average ferromagnetic field experienced by the pair. 
That means, that the spin-flip processes not only modify the oscillation length but also lead 
to an extra decay of the Cooper pair potential and, at a critical strength, these scattering 
processes can completely suppress the 7r-phase superconductivity. 

The scenario of the — 7r-phase transition, as well as the mechanisms of critical current 
enhancement, in the limit of thin F layers differ from those for thick F layers. However, the 
basic physics behind spin-orbit scattering effect for the SF sandwiches where there is no 
induced order parameter oscillations, has not been discussed till now. The purpose of this 
paper is to clarify the mechanisms of influence of spin-orbit scattering on critical current in 
SFIFS tunnel structures with thin F layers. Our analysis is based on the extension of cal- 
culations [6,8,14], so as to include the effects of the spin-orbit scattering. Namely, we study 
the tunnel junction with leads formed by the proximity coupled SF bilayers of a massive 
s-wave superconductor and a thin F metals, when the spin-orbit scattering processes take 
place in ferromagnetic layers. The microscopic model of the proximity effect for SF bilayer 
is employed to discuss the case of arbitrary value of SF boundary transparency, ferromag- 
netic exchange field, and proximity effect strength (Sec. II). Critical current of symmetrical 
SFIFS junction is discussed in Sec. Ill where analytical solutions have been obtained for 
some particular cases. As we shall see, spin-flip scattering plays a major role in trans- 
port properties of superconductor-ferromagnet structures with thin F layers and should be 
considered for understanding the various experimental results. We close with a Conclusion. 



II. SUPERCONDUCTIVITY OF SF BILAYER WITH SPIN-ORBIT 

SCATTERING 



We consider the case when both metals are in the dirty limit condition £$,f > Is,f , 
and the S layer is thick, ds ^> £s, while the F layer is thin, d F <C min(£^, y 'Dp /2nTc). 
Here £s = (D s /2txT c ) 1 ^ 2 and £f = (F> F /2H exc ) 1 / 2 are the coherence lengths of the S and 
F metals; Is and lp are the electron mean free paths, ds,F are the thicknesses, and Ds,f 
are the diffusion coefficients of the S and F metals, respectively. These conditions make it 
possible to neglect the reduction of the critical temperature of the SF bilayer compared to 
that of the bulk S metal, and to imply constant superconducting properties through F layer 
thickness. Throughout this work the F layer will be treated as a single domain film with 
the spin-orbit scattering centers, while there are no such centers in the S layer. Then, the 
superconducting properties of the SF bilayer are described by the Usadel equations [15]. 
The latter can be written as (the domain x > is occupied by the S metal and x < by 
the F metal): 

$ 5CT = A 5 + £ S 2 ^[GWV 5(T ]' , (1) 
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= e^-[Gl^' F± ]' + a so G FT (^ - ^) (2) 

Here Q a = u + iaH exc , a — ± , and = cj„ = nT(2n + 1), n = ±1, ±2, ... is Matsubara 
frequency; aso — 2/3 T so an d tso is the spin-orbit scattering time; the pair potential As 
is determined by the usual self-consistency relations (see, e.g. [16]). We define the x axis 
as perpendicular to the film surfaces and the prime denotes differentiation with respect to 
a coordinate x. The modified Usadel functions [17] = uF S(T /Gsa, &Fa = &aF Fa /G Fa , 
where G Fy s and Fp,s are normal and anomalous Green's functions for the F and S material, 
respectively, are introduced to take into account the normalized confinement of the Green's 
functions G 2 +FF = 1 ; here F(u, H exc ) = F*(cu, —H exc ) , and &sa = ^Fsa/Gsa (see below). 
We also assume that for a nonsuperconducting F metal the bare value of the order parameter 
A F = , however, the Cooper pair correlation function F F ^ due to the proximity effect. 
Eqs. (2) for the F metal are the generalization of the Usadel equations for the case when 
the spin-orbit scattering is present (for more details see Ref. [18]). At temperatures close to 
Tc , when Gs(F)a ~ 1 and $s(F)a/&<r ~ Fs(F)a , Eqs. (2) have been simplified getting the 
form of Eqs. (26) of Ref. [3]. If the spin-orbit processes are absent (tso — ► oo) , the spin 
"up" and "down" subbands do not mix with each other and the equation (2) obtains the 
usual form (see, e.g., [6,14]). 

Eqs. (1), (2) should be supplemented with the boundary conditions. In the bulk of the 
S metal the pair potential is equal to the BCS order parameter A (T) at the temperature 
T : <3>s(oo) = A^oo) = A (T), while at the free (dielectric) interface of the F metal: 
® F ±(— d F ) = 0. Assuming that there are no spin-flip processes at the SF interface, the 
boundary conditions at this interface (see Ref. [14] for details) can easily be generalized for 
the case of two fermionic subband in the form : 

liG F± ^ F± /u ± = isGl^Juj, (3) 



lBFiG F ±<$> F± = u ± G s ±($ s ±/u - $F±Mfc) (4) 

The parameters 7 and ^bf involved in these relations are given by 7 = Ps£,s/pf£, , Ibf = 
Rb/pf^ , where ps,F are the normal state resistances of the S and F metals, Rb is the 
product of SF boundary resistance and its area. In Eqs. (1) and (3), and below we have 
used the effective coherence length £ = (D F /27rT c ) 1 ^ 2 thus providing the regular crossover to 
both limits Tq ^> H exc — > and H exc ^> Tq. The relation (3) provides the continuity of the 
supercurrent flowing through the SF interface at any value of the interface transparency for 
spin "up" and "down" fermionic subband separately, while the condition (4) accounts for the 
quality of the electric contact. An additional physical approximation for the relations (3) and 
(4) to be valid is the assumption that spin discrimination by the interface is unimportant, 
i.e., the interface parameters involved are the same for both spin subband. Generalization 
of Eqs. (3) and (4) to the SF interface with different transmission probabilities for "up" 
and "down" spin quasiparticles is straightforward, however the case we consider contains 
all the new physics we are interested in and is simpler. In such ferromagnetic metals as 
Ni, Gd, etc. the polarization of the electrons at low temperatures is not more than 10%, 
and one can expect the model under consideration reflects the transport properties of these 
ferromagnets hybrid structures. 
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Due to a small thickness of the F metal the proximity effect problem can be reduced to 
the boundary problem for the S layer and a relation for determining <&p± at x = . There 
are three parameters which enter the model: 7m = jdp/£ ? 1b = 7rfgW£ and the energy 
of the exchange field H exc . Using the system of equations (3) and (4), one can obtain the 
equations determining the unknown value of the functions $^ ± (0) and boundary conditions 
for <3>s± . Due to these boundary conditions all the equations for the functions §s± and $f± 
are coupled. In the general case, the problem needs self-consistent numerical calculations 
(similar to those as, e.g. , in Ref. [19]). Here, however, we will not discuss the quantitative 
solution, but pay attention to the qualitative one to consider new physics we are interested 
in. 

We will hereinafter assume that aso <S 1 and dp/^ <C 1 and solve differential equations 
by an iteration procedure finding the corrections to $f±(x) and $>s±{%) m small parameters 
dp/^ and aso- We also restrict ourselves to the quite realistic experimental case of a bilayer 
with a weak proximity effect and low transparency of the SF boundary, i.e., 7m <C 1 and 
7b > 1 . Then the problem is simplified (see Ref. [18] for details) and reduced to the Usadel 
Eqs. (1) for the S layer with boundary conditions (3), while Eq. (4) reduces to 

£,sG s ±%±\x=o ~ 1M^± {1 =F 2a so -^-z ( — —kT c + tj-)} *=o, ( 5 ) 

7tT c A± uj + uj_ uj 7b^± u A± 

where we abbreviated A± — [1 + 2*f B Gsw± + 7|d> 2 _] 1/ ' 2 , G s = 00/ (oo 2 + Aq) 1/2 and jb(m) = 
1b{m)/kTc (functions that are multiplied by aso have been used in the limit aso 0). 
Now, the equation for F layer function $f±(0) has the form: 

$ F± (0) » G 5± $ 5± {1 ± 2« 50 ^ £ }/[^(7i? + G s± /a) ± ]U = o (6) 

In zeroth approximation in 7m the $g ± (0) = 0. So, we can neglect the suppression of 
superconductivity in the S layer assuming that 3>s±(x) is spatially homogeneous: &s±{ x ) = 
A s (x) = A (T) . In the next order in 7 M , by linearizing the Usadel equations for the 
3>s±(x) and making use of the relation (3), the general solution of the linearized equation 
(1) is given by 

* s± (x) = A {1 - C± exp(-/3a;/&)} (7) 

where /3 2 = (c<j 2 + A 2 ,) 1 / 2 /^!^. Substituting this solutions into the boundary relations (5), 
we get for C± : 



C ± - _ 4 {1 T- 2a 50 ^(^ + (8) 



Using Eqs. (7) for x = and relations (6), one can find expressions for $ f± (cj, 0). However 
we will not show here these expressions because of their cumbersome structure. For H exc — > 
the quantity uo± — > , and solution (7), (8) reproduces the results obtained in Ref. [16] for 
an SN bilayer. If tso — » 00, the expressions restore earlier results for SF bilayer (see Refs. 
[6,14]). 

Let us give some comments on the results obtained. For the structure under considera- 
tion, one can expect a kind of induced magnetic properties for S layer, that are the result of 
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the effect similar to the superconducting proximity effect [20]. In accordance with Eqs. (7), 
(8), we see that due to proximity induced magnetic correlations even in the S layer fermionic 
symmetry of subbands has been lost and Cooper pair mixes with its spin-exchanged coun- 
terpartner. 



III. CRITICAL CURRENT OF SFIFS TUNNEL JUNCTION 

We assume that both banks of the Josephson SFIFS tunnel junction are formed by 
equivalent SF bilayers, and the transparency of the insulating layer is small enough to 
neglect the effect of the tunnel current on the superconducting state of the electrons. The 
plane SF boundary can have arbitrary finite transparency, but it is large compared to the 
transparency of the junction barrier. The transverse dimensions of the junction are supposed 
to be much less then the Josephson penetration depth, so all quantities depend only on 
a single co-ordinate x normal to the interface surface of the materials. Using the above 
obtained results we investigate the influence of the spin-flip scattering on critical Josephson 
current in SFIFS tunnel junction. The critical current of the (SF)lI(FS)r tunnel contact 
can be written in the form (see, e.g. Ref. [6]): 

Jc = (eR N /27rT c )I c = —Re > \ L \ R (9) 

where Rn is the resistance of the contact in the normal state; the subscript L (R) labels 
quantities referring to the left (right) bank and the sign of the exchange field depends on 
mutual orientation of the bank magnetizations. 



A. Parallel orientation of the layer's magnetizations 



For parallel alignment of the layer's magnetizations, i.e., with ul = Or , the expression 
for critical current reads: 

ft i o- ri -2-2 , a - iHexc G 2 s <&Sa®S(T -, -l 

{1 + 2-f B GsaVa + -f B UJ a + (TAasolB-^-: ~ } (10) 

We begin with an analytical consideration for the case of a vanishing effective pair breaking 
parameter near the SF boundary, 7m = 0; i.e., the influence of the F layer on the supercon- 
ducting properties of the S metal can be neglected and the order parameter in the S region 
is spatially homogeneous: $s±(x) = A (T) (see, Eqs. (7) and (8)). For the amplitude of the 
Josephson current we than obtained 

fc M - 2^ £ ^^{1 + 2^G S + - Hl c ) + ^sofsj^r} x 

uj>0 exc 

{[1 + 2^ B uGs + lUcu 2 - H 2 exc )f + AH 2 exc f B {G s + ^OT 1 (11) 
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If H exc — > the expression (11) restores the result for SNINS junction (see, e.g., Eq.(28a) 
of Ref. [17] for parameters value under consideration). As is seen from this expression, for 
large enough H exc , the supercurrent changes its sign, i.e., with increasing magnetic energy 
the junction crosses over from 0-phase type to 7r-phase type superconductivity. However, 
the spin-flip processes exert influence upon junction's tendency to set in the 7r-phase state. 

On Fig. 1 we plot a family of the Josephson current amplitude (11) as function of H exc 
when 7^ = 2 , 7m = for various values of the spin-orbit scattering intensity aso ■ Here we 
also show the function 5j^ M = jc M ( a so) ~ jc M ( a so — 0) .The main feature here is that 
the increase of the exchange energy pulls the SFIFS junction to 7r-state. The new result 
of this figure is that, as intensity of spin-orbit scattering processes increases, the critical 
current amplitude, jc< M ( a so), for the 7r-state decreases. The spin-orbit scattering reduces 
the effect of exchange field, and this reduction has nonlinear character with its maximum in 
the region where transition to the 7r-state take place. 

A more realistic case is shown on Figs. 2, where we plot the results of numerical cal- 
culations of the critical current, Eq. (10), for the tunnel junction when suppression of the 
S layer superconductivity occurs due to proximity effect (7m ^ 0). The function 5j^ M has 
also been shown. Again, the main result of these calculations is that the spin-flip processes 
can sizably reduce the SFIFS junction tendency to 7r-phase state. In contrast to a simple 
physical picture, the suppression has nonlinear behavior getting its maximum near H exc 
region of cross over from 0-type to 7r-type superconductivity. While our results have been 
obtained for the SF bilayers with a weak proximity effect and low boundary transparency, 
qualitatively the conclusions are valid for arbitrary 7# and 7m values. 



B. Antiparallel orientation of the layer's magnetizations 

To be definite, we took Cjl = oj + iH exci ujr = u — iH exc . Now the critical current of the 
(SF) l I(FS)r tunnel contact can be given in the form 

•AF o T p Q V^ G S+ $s + G S -$S- un ,- ~ n2 , Gl+^S+^S+f . iHexc ^ 

3c =2 T^ Re 2-. -2 {(G s+ + 1bco + ) + (i + 4a 50 ^-)]} 

x{(G S - + 7^-) 2 + ggzg*pggz (i _ Aaso^)}- 1 (12) 

Fist we assume that one can neglect the suppression of the superconductivity in the S layer 
(7m = 0). Then, after simple but cumbersome algebra we have for the critical current 

Wl c f B (Gs + lBU?r in (13) 

For H exc — > expression (13) restores the result for SNINS junction (Eq. (28a) [17]). One 
can see, that for uj < H exc the expression in figured brackets is low than one. As a result, the 
0-phase state across the contact has been preserved, and in some interval of the exchange 
field energy the enhancement of dc Josephson current occurs. The spin-orbit processes exert 
upon this behavior, reducing the current enhancement. 
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On Fig. 3 we show a family of the Josephson current amplitude (13) as function of 
H exc when 7# = 2 , 7^ = for various values of the spin-orbit scattering intensity aso ■ 
Here we also plot the function 8jq F = jc F { a so) ~ jc F ( a so — 0) .The main and unusual 
feature here is that the increase of the exchange energy enhances the SFIFS junction 
critical current. The new result of this figure is that, as the intensity of spin-orbit scattering 
processes increases, the tendency of the jc F ( a so) to enhance decreases. The spin-orbit 
scattering reduces the effect of exchange field, and this reduction is nonlinear in character 
with its maximum in the most interesting for experimentalists region. More general cases 
are shown on Fig. 4, where we illustrate the results of numerical calculations of the critical 
current, Eq. (12), taking into account a suppression of the S layer superconductivity by 
the proximity effect ('Jm 7^ 0). The function Sj^ F has also been shown. One can see again, 
that the spin-orbit suppression has nonlinear behavior getting its maximum in the most 
interesting for experimental investigation region of H exc . Qualitatively, these results are 
valid for arbitrary values of 7^ and 7m • 

IV. CONCLUSION 

We do not consider here the experimental situation, which is now unclear an controver- 
sial even for the SF hybrid structures with thick ferromagnetic layers. Let us only note, 
that spin-orbit scattering is relevant for ferromagnetic conductors containing large Z num- 
ber. Magnetic inhomogeneity of the materials, such as the F layer multidomain structure, 
domain walls, inhomogeneous " cryptoferromagnetic" state imposed by superconductor also 
give nonzero probability amplitude for the spin-flip scattering. For nanoscale hybrid struc- 
tures strong electric field arising near metal-metal boundaries is also an important source 
of the spin-flip processes. In the absence of any precise information about the magnetic 
structure of the samples used in experiments on SF sandwiches, we restrict ourselves to 
making only qualitative calculation. 

In conclusion, we investigate the spin-orbit scattering effect on critical current of SFIFS 
tunnel junction for the case of thin F layers, when the superconducting order parameter 
oscillation is absent. Instead, the parameter phase jumps at the SF interfaces. The analysis 
is based on a microscopic theory for proximity coupled S and F layers. The main result 
of our calculations is that the spin-flip processes can sizably modify the behavior of the 
dc Josephson current versus the exchange field for SFIFS tunnel junction. We find that 
the reduction of the exchange field effects is nonlinear in character getting its maximum 
in the ferromagnetic field's region where the critical current enhancement or the transition 
to the 7r-state occur. Hence, for understanding the various experimental results on tunnel 
structures with thin F layers, the coupled effects of the exchange interaction and spin-orbit 
scattering must be considered. 

The authors are grateful to Valery Ryazanov and Elena Koshina for useful discussions. 
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Figure Captures 

Fig. 1. Critical current of SFIFS tunnel junction with parallel orientation of the F 
layers magnetization versus exchange energy for various aso/^o — 0, 0.05 , 0.1 and 
0.15 ; the SF interface transparency is low 7# = 2 and proximity effect is weak 7^ = 0. The 
additive part of critical current due to spin-orbit scattering, Sj^ M , is also shown. T = O.lTc 

Fig. 2. Critical current of SFIFS tunnel junction with parallel orientation of the F 
layers magnetization j^ M versus exchange energy for various 7m = 0, 0.05, 0.1, 0.15; 7^ = 2 
and aso/^o — 0.1. The additive part of critical current due to spin-orbit scattering, 8jq 
, is also shown. T = O.IT^ . 
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Fig. 3. Same as on Fig.l but for antiparallel orientation of the F layers magnetization. 
The additive part of critical current due to spin-orbit scattering, Sj^ F , is also shown. 

Fig. 4. Same as on Fig. 2 but for parallel orientation of the F layers magnetization. The 
additive part of critical current due to spin-orbit scattering, 5j£ M , is also shown. 



10 




Fig. 2 




Fig. 3 




Fig.4 



